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The electron-phonon coupling in potassium-doped graphene on Ir(lll) is studied via the renor- 
malization of the 7r* band near the Fermi level, using angle-resolved photoemission spectroscopy. 
The renormalization is found to be fairly weak and almost isotropic, with a mass enhancement 
parameter of A = 0.28(6) for both the K — M and the K — Y direction. These results are found to 
agree well with recent first principles calculations. 



Many-body effects in graphene have recently attracted 
considerable attention, mainly because of graphene's per- 
ceived role as a future electronics material, but also be- 
cause graphene is highly suited as a fundamental test case 
for the investigation of many body effectsi*2a2aiaSaSaLS4&. 
In particular, the structural and electronic simplicity 
of graphene seems to imply that it is accessible to 
both first-principles calculations and angle-resolved pho- 
toemission (ARPES). Other advantages are that high- 
quality graphene can be prepared and doped by an elec- 
tric field or by adsorption. 

This paper is concerned with the electron-phonon in- 
teraction in electron-doped graphene. Despite the above- 
mentioned apparent advantages of graphene to study 
many-body interactions, there has been considerable dis- 
pute on the strength of the electron-phonon coupling in 
both experimental and theoretical studies. From the ex- 
perimental side, there are several possible causes for this. 
First of all, experiments have not been performed on 
free-standing graphene and it is not easy to evaluate the 
role of the substrate. More importantly, determining the 
electron-phonon self-energy from the dispersion's renor- 
malization is not simple, after all, since the unrenormal- 
ized (bare) dispersion is also influenced by many-body 
effects at higher energies^. The usual approach of us- 
ing a linear bare dispersion fails for graphene 3 , a fact 
which has lead to unrealistic results in early studies 4 . 
Temperature-dependent studies, usually a stable alter- 
native to determine the coupling strength in the form of 
the electron-phonon mass enhancement parameter A, are 
of limited value because of graphene's exceptionally high 
Debye temperature. At this point, rather different ex- 
perimental results for the size of A and its variation over 
the Fermi surface have been reported for graphene and 



graphite intercalation compounds 4 ^ 1Q i 11 . 

Theoretical results generally point towards a rather 
weak coupling but the absolute values of A vary greatly 
between different calculations 3 -^. Moreover, it has been 
pointed out that the local density approximation might 
not be able to predict the electron-phonon coupling 
correctly^. On the positive side, a very recent calcula- 
tion of the expected linewidth in ARPES, including both 
electron-electron and electron-phonon scattering effects 
gives good agreement with experimental data for doped 
graphene on SiC, at least in the K — M direction for 
states near the Fermi energy (Ef) 6 . 

In this paper, we analyze the electron-phonon coupling 
strength for electron-doped graphene grown on Ir(lll). 
This approach has the potential advantage of starting 
with a graphene layer which is relatively well decoupled 
from the substrate and only very weakly doped when 
clean. We use a stable method to determine A despite 
the unknown bare dispersion. We find a moderate A 
value which varies little over the Fermi surface, a scenario 
that agrees well with a recent theoretical prediction 3 . We 
compare our results to other works using graphene on SiC 
as a starting materiali^ as well as the KCg and CaC6 in- 
tercalation compound ^ 10 ! 11 . In KCs, and probably also 
in CaC6, the graphene sheets are effectively decoupled 
because of the increased distance between them after al- 
kali intercalation 13 . 

Experiments were performed at the SGM3 beamline 
of the ASTRID synchrotron radiation facility at Arhus 
University^. The endstation is equipped with a 150 mm 
hemispherical electron energy analyzer (Specs). The data 
were taken with a photon energy of 45 eV, and at a sam- 
ple temperature of 70 K. The total energy and k resolu- 
tion amounted to 28 meV and 0.01 A -1 , respectively. 
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The graphene film was prepared on Ir(lll) using a 
standard recip o 15 i 16 . The cleanness of the Ir substrate be- 
fore graphene formation was monitored using the 4/ sur- 
face core level shift, and the quality of the graphene film 
was checked by means of Low Energy Electron Diffraction 
(LEED) and ARPES, which showed the characteristic 
mini-gaps close to the K points. Potassium was evapo- 
rated from a commercial getter source (SAES). The data 
reported here were taken in a situation in which LEED 
showed a (2 x 2) superstructure. This corresponds to 
the KCg phase in an alkali metal / graphite intercala- 
tion compound. At this doping level, the Dirac point is 
shifted to a binding energy of 1.29 eV and the electron 
concentration is ~ 1 x 10 14 cm -2 or, equivalently, the 
doping level is 0.054 extra electrons per graphene unit 
cell. 

An area of two-dimensional reciprocal space around 
K was sampled and a resulting dispersion and the Fermi 
contour are shown in Fig. [T](a) and (b), respectively. De- 
tailed dispersions along the K — M and K — Y directions 
were extracted and shown in Fig. [TJc) and (d), respec- 
tively. The dispersion along the K — Y direction in (c) is 
a sum of four parallel, closely spaced slices. Along K — M 
only a single cut has been used. Averaging over several 
parallel cuts in this direction gives rise to broadening, as 
expected because of the sharp corner of the Fermi surface 
in this direction. The dispersions along K — M and K — Y 
both show a pronounced renormalization of the bands, as 
evident from the kink in the dispersions around a binding 
energy of « 180 meV. The apparent strength of the kink 
is greater in the K — M direction than in the K — Y direc- 
tion. Qualitatively, these results are in good agreement 
with previous findings in similar systems 4 ^ 10 ! 11 . 

Due to the three- fold symmetry of the Fermi surface, 
several equivalent cuts can be made. For the following 
analysis we have focused on the directions with the high- 
est signal to noise ratio and smallest contribution of the 
underlying Ir surface states. The results from other direc- 
tions are consistent with the reported findings and have 
been used to estimate the uncertainty of the obtained 
coupling strengths. 

The detailed renormalized dispersions were extracted 
by fitting Lorentzian lines with a linear background to 
the momentum distribution curves (MDCs), i.e. to con- 
stant energy profiles. The fits were of good quality even 
though, strictly spoken, MDC cuts through the spectral 
function can only be described by Lorentzian lines in the 
case of a linear dispersion. This procedure also gives the 
energy-dependent full width at half maximum (FWHM) 
of the MDC peaks. 

In the next step of the analysis, the self-energy E has 
to be extracted from the data. E is a complex quan- 
tity and its real and imaginary parts are related via a 
Kramers-Kronig transformation. The measured spectral 
function contains information about both the real part E' 
and the imaginary part E" of the self-energy via the size 
of the renormalization and the linewidth of the MDCs, 
respectively^. Extracting E" from the MDC width suf- 



r < — k — ► m 




1-5 J 1 ^ , r- 

-0.2 0.0 0.2 0.4 
k(A- 1 ) 




-0.2 -0.1 0.0 0.1 -0.1 0.0 0.1 0.2 



k-k F (A- 1 ) 

FIG. 1: Photoemission intensity around the K point shown as 
different cuts through a three-dimensional data set. (a) large- 
scale dispersion through K. (b) Fermi surface around K with 
the dashed line indicating the cut shown in (a). More detailed 
Fermi energy crossings along the K — Y (c) and K — M (d) 
directions, (e) and (f) show the simulated photoemission in- 
tensity corresponding to the data in (c) and (d), respectively. 



fers from the problem that the width is usually rather 
noisy and that, strictly, E" can only be determined in 
the case of a linear bare dispersion. Using the renor- 
malized dispersion to extract E', on the other hand, 
requires the knowledge of the bare dispersion. In the 
present case, both restrictions pose formidable problems 
but several approaches have been suggested in order to 
overcome them, such as measuring the dispersion at high 
temperatures^ or using a self-consistent procedure to ex- 
tract We employ a procedure which is similar to the 
latter approach. 

In the first step of this analysis, a model for the mea- 
sured spectral function is calculated from initial guesses 
for the bare dispersion and E. The bare dispersion is as- 
sumed to be a second order polynomial. E is calculated 
from an Eliashberg coupling function a 2 F(uu), consisting 
of the contributions of 5 Einstein oscillators with ener- 
gies that are evenly distributed over the energy range 
from 21 meV to o; max = 190 meV. The oscillators, with 
energy uoi and a coupling strength value , have a width 
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of 1 meV. From a 2 F(uj), it is straightforward to calculate 
E" (and thus also £') through 



a 2 F(wO[l-/(w-u/)+/(wW)+2n(w')]dw' 

(1) 



where n and / are Bose and Fermi functions, respectively. 
Note that the simple, discrete character of the Eliashberg 
function has only a minor effect on Indeed, reason- 
able fits of E" can sometimes be achieved with even sim- 
pler models^. Adding further Einstein oscillators in the 
present case would render them redundant in view of the 
experimental energy resolution. 

We add an offset to E" in order to account for energy- 
independent defect scattering. As we are only concerned 
with a small energy window close to Ep, we ignore the ef- 
fect of electron-electron scattering altogether. The model 
thus contains 9 variable parameters: 3 governing the 
bare-particle dispersion, 5 values of A^ and one offset for 
Using these, the spectral function can be calculated 
from the bare dispersion and E. The result is multiplied 
with a Fermi function and convoluted with the experi- 
mental energy and k resolution functions. 

In a second step, the MDC position and width are ex- 
tracted from this simulated spectral function in the same 
way as from the experiment. The MDC peak positions 
and widths of experiment and simulation are compared 
using a combined \ 2 f° r position and width. The model 
parameters are then optimized using a steepest-descent 
algorithm. 

The results of this fitting procedure are given in Fig. [2] 
and[3j Fig. [2]shows the agreement between data and sim- 
ulation as well as the bare dispersion obtained from the 
fit. Fig. Ogives the resulting E. The level of agreement 
in Fig. [2] is clearly satisfactory given the simplicity of 
the model. The self-energy obtained from the procedure 
reveals the expected contributions caused primarily by 
the high energy optical phonon modes and, less strongly, 
lower-lying modes, consistent with the literature ^ 1Q i 11 . 
The mass enhancement parameter A can now be calcu- 
lated as the first reciprocal moment of a 2 F{uo). We find 
the same A value of 0.28(6) for the K - M and K - T 
directions. 

It is quite remarkable that the obtained A values are 
so similar for the two directions despite of the apparently 
stronger renormalization along K — M. The reason is the 
downward curvature of the band in this direction as it ap- 
proaches the van Hove singularity, as noted earlier—. This 
curvature also leads to counter-intuitive slight increase of 
the linewidth upon approaching the Fermi level, best ob- 
served between 450 and 200 meV in Fig. EJd) and also re- 
ported for similar compounds 10 . An energy dependence 
of the linewidth in regions with a rapidly changing den- 
sity of states is, of course, expected and has been found 
for other systems^ -. The similarity of the A values can 
already be guessed from inspecting E" alone. The abso- 
lute value of the two curves in Fig. [3] is different but the 
increase in E" due to electron-phonon coupling is very 
similar. Note that there is no a priori reason to assume 
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FIG. 2: (Color online) Detailed comparison between the mea- 
sured spectral functions and those found in the iterative op- 
timization procedure outlined here, (a) and (c) peak position 
in MDC cuts together with the bare dispersion obtained from 
the fit (dashed line); (b) and (d) peak full width at half max- 
imum. Points with error bars represent data, solid lines the 
fit. 



anisotropic A values. If the coupling to very low-energy 
acoustic phonons (i.e., with q vectors significantly shorter 
than the size of the K Fermi contour) were relevant, one 
could expect shorter lifetimes around the K — M direc- 
tion, as more states with similar energy are close by. But 
as coupling to high energy optical phonons is obviously 
dominant and intervalley scattering is thought to be im- 
portant for electron-doped graphen o 1Q i n , this simple ar- 
gument does not hold. 

We note that the outlined approach for determining E 
and A has several advantages over other methods. First 
of all, it does not rely on any knowledge of the bare dis- 
persion, as long as it is locally parabolic. In the case of 
graphene, this is crucial, because the strong non-linearity 
in the bare dispersion close to the van Hove singularity 
does not only upset a simple determination of E', it also 
makes it impossible to relate the MDC linewidth to E" 
through the usual relation that E" = vFWHM/2, where 
v is the group velocity of the state and FWHM the MDC 
full width at half-maximum. As mentioned before, this 
effect is most pronounced along K — M. 

The data analysis approach also has the advantage of 
built-in consistency between E' and E" at the expense of 
assuming particle-hole symmetry for the transformation 
from £" to £'. Also, it permits sufficient variation of E, 
while describing the system with a manageable number 
of parameters. Indeed, the idea of using an Eliashberg 
function constructed with a small number of Einstein os- 
cillators is simple and stable and one might wonder why 
it works so well. The reason is that one does not fit the 
actual Eliashberg function but E, or rather the observed 
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FIG. 3: (Color online) Optimized self-energies (real part X' 
and imaginary part X") for the K — Y (a) and K — M (b) 
directions. 

dispersion, which contains only a temperature-broadened 
integral over the a 2 F(uj). Details in the latter become in- 
significant but the determination of the coupling strength 
is stable. Finally, the procedure includes the band distor- 
tion caused by the finite resolution of the spectrometer. 
In particularly the energy resolution is well-known to dis- 
tort the dispersion close to Ep, leading to an error in the 
band slope 18 . For the determination of A this is obviously 
very significant because for T = 0, A is equal to the slope 
ofE'atEjr. 

We compare our results to A values obtained for 
graphene and graphite with similar electron doping. For 
KCs it has been recently shown that the very interca- 
lation of K atoms also leads to an effective reduction of 
the system's dimensions from three to two, caused by the 
larger distance between the graphene sheets, such that a 
direct comparison to graphene is possible 11 . 

The analysis of the electron-phonon coupling for 
electron-doped graphene grown epitaxially on SiC has 
given rather similar A values as reported here 1,2,4 . An 
early study had reported a rather large anisotropy, espe- 
cially at even higher doping levels, which was most proba- 
bly caused by the unrealistic assumption of a linear bare 
dispersion^. On KCs, on the other hand, A was found 
to be rather anisotropic and also stronger than here 11 . 
The reason for this is unclear and could lie in the differ- 
ent approach to data analysis or in the fact that having 
K atoms in between graphene sheets leads to markedly 
different phonon dispersion, even though the electronic 
properties of the individual graphene sheets are similar to 
free graphene. Our results agree well with a recent first 
principles calculations of the electron-phonon coupling 



in electron-doped graphene^ which predicts an almost 
isotropic electron-phonon coupling of similar strength. 

Finally we compare the absolute MDC linewidth close 
to the Fermi level to that in similar systems in order to 
assess the role of defect scattering in this system. The 
linewidth found here is roughly 0.095 A -1 which is quite 
similar to data from KCg 11 but significantly broader than 
for alkali-doped graphene on SiC 1 . It has to be noted, 
however, that the graphene/ SiC data in Ref. 1 have been 
taken at lower doping levels and MDCs taken for the 
highest coverages suggest that there is an increase of 
linewidth in this regime. A plausible explanation for an 
increased linewidth for graphene on Ir(lll) is the fact 
that the resulting superstructure is incommensurate^!, 
implying a loss of translational symmetry for the com- 
bined system and a presumably increased scattering rate, 
even if the coupling between substrate and graphene is 
small. We also note that the adsorption of alkali metals 
on graphene/ SiC and graphene/Ir(lll) appears to be dif- 
ferent. On SiC it is possible to continuously vary the dop- 
ing level whereas this cannot be done for Ir(lll). Here 
the potassium adsorbates assemble in islands with (2 x 2) 
periodicity and at lower coverages two tt bands are ob- 
served, one similar to the doped case reported here and 
one similar to clean graphene on Ir(l 11)22. Finally, it 
is questionable if a comparison of the coupling strengths 
and linewidths for different doping levels is at all mean- 
ingful. Valla et al. have evoked a dynamical inter valley 
nesting effect to explain the observed coupling on CaC6 
and this effect implies a strong doping dependence of the 
coupling 10 . 

In conclusion, we have determined the electron-phonon 
coupling strength for a layer of alkali-doped graphene on 
Ir(lll) in the K — V and K — M directions. We find 
an almost isotropic scattering strength consistent with 
recent first-principles calculations. The approach used 
here takes the non-linearity of the bare dispersion and 
the finite experimental resolution into account and uses 
a simple model for the Eliashberg function which allows 
a stable, reliable and self-consistent analysis of the mea- 
sured spectral function. 
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